Anisotropy induced Feshbach resonances in a quantum dipolar gas of magnetic atoms 
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We explore the anisotropic nature of Feshbach resonances in the coUision between ultracold mag- 
netic submerged-sheU dysprosium atoms, which can only occur due to couplings to rotating bound 
states. This is in contrast to well-studied alkali-metal atom collisions, where most Feshbach reso- 
nances are hyperfine induced and due to rotation-less bound states. Our novel first-principle coupled- 
channel calculation of the collisions between open-4f-shell spin-polarized bosonic dysprosium reveals 
a striking correlation between the anisotropy due to magnetic dipole-dipole and electrostatic in- 
teractions and the Feshbach spectrum as a function of an external magnetic field. Over a 20 mT 
magnetic field range we predict about a dozen Feshbach resonances and show that the resonance 
locations are exquisitely sensitive to the dysprosium isotope. 
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A strongly interacting quantum gas of magnetic atoms, 
placed in an optical lattice, provides the opportunity 
to examine strongly correlated matter, creating a plat- 
form to explore exotic many-body phases known in solids, 
quantum ferrofluids, quantum liquid crystals, and super- 
solids PJI24. Recent experimental advances [BHTO] in trap- 
ping and cooling magnetic atoms pave the way towards 
these goals. 

In general, interactions between magnetic atoms are 
orientationally dependent or anisotropic. At room tem- 
perature anisotropic interactions are much smaller than 
kinetic energies and other major interactions between 
atoms, therefore can be ignored. The situation is differ- 
ent for an ultracold gas of atoms with a large magnetic 
moment. It was, for example, demonstrated that the 
anisotropy due to magnetic dipole-dipole interactions be- 
tween ultracold chromium atoms leads to an anisotropic 
deformation of a Bose Einstein condensate (BEC) [TT]. 
Moreover, anisotropy plays a dominant role in coUisional 
relaxation of ultracold atoms with large magnetic mo- 
ments [SHZIIIMI]- 

In this Letter we pursue ideas for using anisotropic 
magnetic and dispersion interactions to control collisions 
of ultracold magnetic atoms by using Feshbach reso- 
nances [in]. Resonances, shown schematically in Fig. [l] 
appear when the energy of "embedded" bound states 
cross the energy of the entrance channel or initial scat- 
tering state. The embedded state is a level of a potential 
dissociating to a closed channel whose asymptotic energy 
is larger than that of the entrance channel. Coupling with 
the entrance channel leads to a resonance. 

Feshbach resonances make it possible to convert a 
weakly interacting gas of atoms into one that is strongly 
interacting and along the way promise to make avail- 
able many of the collective many-body states mentioned 
above. Alternatively, interactions can be turned off all 
together to create an ideal Fermi or Bose gas, for which 



thermodynamic properties are known analytically. Fes- 
hbach resonances can also be used to create BECs of 
weakly-bound molecules pTfj , which can be optically sta- 
bilized to deeply-bound molecules [TH]- For fermionic 
atoms the BCS-BEC phase transition [12] and univer- 
sal many-body behavior of strongly interacting magnetic 
atoms can be studied via Feshbach resonances. Finally, 
three-body Efimov physics [20] can be explored. 
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FIG. 1; (Color online) Potential energy curves for a 
^®*Dy-|-^^*Dy collision in an external magnetic field B as a 
function of internuclear separation. The (red) dashed line 
with zero energy indicates the energy of the entrance chan- 
nel. Two Feshbach resonances are schematically shown by 
(red) horizontal lines, which end at the classical outer turn- 
ing point of a closed channel. Their energy increases, in- 
dicated by arrows, with magnetic field and a resonance oc- 
curs when this energy equals the entrance-channel energy. 
There are 91 diagonal potential matrix elements for chan- 
nels \{jiJ2)jmj,£me} with nij + mi — —16 and even £ < 10. 
We use B = 50 G. The curves are colored by their m,j value, 
while for mj = — 16 curves their £ value is indicated. Here 
1 G=0.1 mT, an = 0.0529177 nm is the Bohr radius, and 



k = 1.38065 ■ lO"^"* J/K is the Boltzmann constant. 



The most promising atoms to look for the effect of 
anisotropy on coUisions are submerged-sheU atoms, which 
have an electronic configuration with an unfilled inner 
shell shielded by a closed outer shell. In particular, we are 
interested in the rare-earth dysprosium (Dy) atom with 
a ^Ig ground-state, a total angular momentum j — 8, and 
a large magnetic moment of ~ lO/is, for which the elec- 
tron spins of the inner 4f^° shell are aligned such that 
its orbital angular momentum is maximal and largely 
unquenched. Here ^b is the Bohr magneton. As a 
result, Dy's magnetic and electrostatic properties are 
highly anisotropic. A quantitative description of the col- 
lision between two dysprosium atoms is challenging. For 
example, our previous study [15] showed that there are 
153 Born-Oppenheimer potentials that dissociate to the 
ground ^Ig -|- ^Is state. 

We present a first-principle coupled-channel model 
allowing us to calculate anisotropy-induced magnetic 
Feshbach-resonance spectra of bosonic Dy atoms. The 
model treats the Zeeman, magnetic dipole-dipole, and 
electrostatic isotropic and anisotropic dispersion inter- 
actions equally. Bosonic Dy isotopes have zero nuclear 
spin. Thus, there is no nuclear hyperfine structure and 
only Zeeman splittings remain. The weak quadrupole- 
quadrupole interaction [15J is included for completeness. 

We already note that Feshbach resonances in rare- 
earth magnetic atoms are different in nature than those 
in alkali-metal atom collisions. As pointed out in Ref. [15] 
for a coupled-channel calculation there exists at most one 
channel with zero relative nuclear orbital angular mo- 
mentum £, which for ultra-cold collisions is also the en- 
trance channel. Consequently, resonances occur due to 
anisotropic coupling to bound states with non-zero £. 

The focus of this Letter is on ultra-cold collisions of 
atoms prepared in the energetically-lowest Zeeman state 
j = 8 and projection m = —8. Inelastic exothermic 
atom-atom processes, where the spin projection of one or 
both of the atoms changes, are absent and, consequently, 
Feshbach resonances can be readily observed. 

We start by setting up the Hamiltonian, interatomic 
potentials, and channel basis for two bosonic ^Ig Dy 
atoms with zero nuclear spin. This Hamiltonian assum- 
ing a magnetic field B along the z direction is 

where R describes the orientation of and separation be- 
tween the two atoms. The first two terms are the ra- 
dial kinetic and rotational energy operators, respectively. 
The Zeeman interaction is Hz — gsfJ'BiJiz + J2z)B with 
gs — 1.24159 the g-factor of Dy [3T] and jiz is the z com- 
ponent of the angular momentum operator ji of atom 
i = 1,2. The electronic Hamiltonian, including nuclear 
repulsion, V{R, r) is anisotropic and r labels the elec- 
tronic variables. Finally, jjr is the reduced mass and for 
i? — >■ oo the interaction V{R, r) — )■ 0. 



Our coupled-channels calculations [22] are per- 
formed in the atomic basis |(jij2)jWj,^mf) = 
YimiiO A)\{iiJ2)j'mj) , where Yfmf (6*, 0) is a spherical har- 
monic and angles 9 and (j) give the orientation of the in- 
ternuclear axis relative to the magnetic field direction. 
In this basis the Zeeman and rotational interaction are 
diagonal with energies gsfJ.B'mjB + H^£{£ + l)/{2^rR^)- 
Coupling between the basis states is due to V{R,t) and 
wiU be discussed in detail below. Excited atomic states, 
for example those with ji ^ 8, are not included as their 
internal energy is sufficiently high that the effects of cou- 
pling to these states is negligible. The Hamiltonian H 
conserves Mtot ~ rrij + mg and is invariant under the 
parity operation so that only even (odd) £ are coupled. 
For homonuclear collisions only basis states with even 
j + £ exist. Figure [l] shows an example of the long-range 
diagonal matrix elements in the atomic basis of the sum 
of the rotational, Zeeman, and electronic Hamiltonian. 
We have used Mtot — —16 and even £ < 10. In fact, 
only the potentials dissociating to the six energetically- 
lowest Zeeman states are shown. The large number of 
potentials indicates the large number of resonances that, 
in principle, are possible. 

Coupling between basis states is due to V{R,t). It is 
convenient to first evaluate this operator in a molecular 
basis with body-fixed projection quantum numbers de- 
fined with respect to the internuclear axis. We use the 
molecular basis \{jiJ2)j^) with projection O of J along 
the internuclear axis. For Dy2 the matrix elements of 
V{R,t) conserve the projection fl but not the length 
j. The eigenenergies of V{R,t) at each value of R are 
the adiabatic (relativistic) Born-Oppenheimer potentials 
[531 mi- Typically, these potentials Un\n\aiR) are ob- 
tained from an electronic structure calculation and la- 
beled by ?t.|^|ct I where |r2| is the absolute value of SI, 
a — g/u is the gerade/ungerade symmetry of the elec- 
tronic wavefunction, and n — 1,2,- ■ ■ labels curves of the 
same \fl\ , in order of increasing energy. For bosonic 
Dy2 the 81 gerade states are superpositions of even j, 
while the 72 ungerade states are superpositions of odd j. 

For R > 27ao, beyond the Le Roy radius where the 
atomic electron clouds have negligible overlap, we assume 
that V{R, t) is the sum of the magnetic dipole-dipole, 
Vtifj.{R) oc 1/-R"^, the electric quadrupole-quadrupole, 
Vqq{R) oc 1/R^, and the van-der-Waals dispersion 
Vdisp{R) oc 1/R^ interaction. Reference [15] reported the 
matrix elements of the operator Vdisp (R) in the molecular 
basis and tabulated the adiabatic Cg.nno- dispersion co- 
efficients obtained by diagonalizing Vdisp (-R). Crucially, 
the cigcnfunctions of Vdisp (^) are independent of R. 

At shorter range coupling between basis states is 
more complex. Rather than determining all Born- 
Oppenheimer potentials, we have opted for the fol- 
lowing approach. First, we calculate the single ger- 
ade potential Uieg{R) with the maximal projection 
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FIG. 2: Elastic rate coefficient of m = —8 ^^*Dy collisions as 
a function of magnetic field using a collision energy of E/k = 
30 nK. Partial waves £ up to 10 are included. 



f2 = 16 using a coupled-cluster method with single, 
double, and perturbative triple excitations (CCSD(T)) 
[25] together with the scalar relativistic Stuttgart 
ECP28MWB pseudopotential and associated atomic 
bases sets (14s,13p,10d,8f,6g)/[10s8p5d4f3g]. The po- 
tential has a minimum at Re = 8.771ao with depth 
Del {he) = 785.7 cm-i. For ^^-^Dya it has a uJe/{hc) = 
25.6 cm~^ and contains 71 bound states. (We omit 
the n — 1 label in Uieg{R)-) We then assume that 
the R < 27ao electronic wavefunctions of the Born- 
Oppenheimer potentials are the same as those deter- 
mined by the dispersion interaction, and the relation be- 
tween energies of the ab-initio potentials is the same as 
for its Cg coefficient. Hence, the adiabatic potentials sat- 
isfy Unna{R)/Un'na{R) = Ca^nQa /Ce,n'n' a for R < 27ao 
and with eigenfunctions as determined by the dispersion 
interaction. Equivalenty, this allows us to write V{R,t) 
in terms of Uieg{R) as 



ViR,T) = V,^{R) + VQQiR) + 



'CQ^lQg/R^ 



-V^ispiR) (2) 



for any R. This approach reduces the number of inde- 
pendent short-range potentials to one. 

For practical reasons it is advantageous to write 
Vdisp(-R) as a sum of spherical tensor operators in the 



TABLE I: Dispersion coefficients c^ in units of Ehaf), where 
Eh = 4.35974 x 10^^** J is the Hartree. The strength of the 
magnetic dipole-dipole and quadrupole-qudrupole interaction 



arec^,^ = -5.0269-10 
respectively. 
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laboratory frame. That is 
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where CkqiO, </>) = ^y■i^^/{2k + l)Ykq{9, (f)) and the spher- 
ical tensors TI^ of rank k and with components q are 



defined by 
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where / is the identity operator and [j ® i']kq denotes a 
tensor product of angular momentum operators J* and / 
coupled to an operator of rank k and component q [26j . 
The higher-order tensor operators are constructed in an 
analogous manner. Equation [4] has three, three, and one 
tensors T,.*' of rank fc = 0, 2, and 4, respectively. The dis- 



persion coefficients 



.(«) 



listed in Table 



strength of each term. Their order in Eq. 
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determine the 
is by decreas- 



ing absolute value. The isotropic Tqq term is the largest 
by far with the strongest anisotropic contribution from 
a dipolar (rank-2) operator constructed from the angu- 
lar momentum of only one atom coupled to the rotation 
of the molecule. Finally, the magnetic dipole-dipole and 
quadrupole-quadrupole interaction are 

Vqq{R) - -^cqqE(-1)«C4.-,(0,0)t4' 



(5) 



(6) 



respectively. Their strengths are listed in Table |T] 

Figure [2] shows the elastic rate coefficient for the col- 
lision between two m = —8 ^^*Dy atoms at a collision 
energy of E/k = 30 nK. We use Mtot = — 16 and include 
channels with even i up to ten leading to a close-coupling 
calculation with 91 channels. Fields up to i? = 200 G fall 
comfortably within the experimentally accessible values 
(16j . The graph shows a dozen of Feshbach resonances; 
some are broad, many are very narrow. By performing 
calculations that include fewer partial waves we have ob- 
served that the resonances can not be labeled by a single 
partial wave. For example, the three broad resonances at 
B « 30 G, 110 G, and 170 G are already present when 
only £ = 0, 2, and 4 channels are included. Their loca- 
tions, however, shift significantly when higher (. channels 
are included and only converge to within a few Gauss 
when £ = 8 channels are included. In general, we find 
that the magnetic-field location of a resonance that ap- 
pears when channels with partial wave £ are included 
stabilizes when channels up to € -I- 4 are included. 



600 



300 




1000 



500 I^J 




600 

300 







J 







50 100 150 

Magnetic field (G) 



200 



FIG. 3: Scattering length of m = —8 ^''^Dy atoms as a func- 
tion of magnetic field with and without the magnetic dipole- 
dipole or the anisotropic contribution of the dispersion inter- 
action. The top panel shows the case when all interactions are 
included. At i? = the scattering length is 89ao. The mid- 
dle and bottom panels are obtained when the dispersion and 
magnetic dipole-dipole anisotropy is set to zero, respectively. 
Even waves £ up to 10 are included. 



We stress that this behavior with increasing number 
of channels is unlike that observed in alkali-metal atom 
collisions [TB] or even in collisions of strongly magnetic 
chromium atoms [27j . For these atoms resonances do not 
shift by more than a few Gauss when additional partial 
waves are added. Hence, resonances can be labeled by 
a partial wave quantum number. For dysprosium the 
anisotropic interactions are so strong that states with 
different partial waves are strongly mixed. We have also 
studied the effect of the uncertainty in the depth of the 
f2 = 16 Born-Oppenheimer potential. The depth De was 
changed by adding a localized correction to UiQg{R) that 
does not modify its long-range potential. A depth change 
by no more than 10 cm~^ changes its number of bound 
states by one. Changing the depth by smaller amounts 
changes the resonance spectrum non-trivially. For exam- 
ple resonance widths can be modified dramatically and 
rate coefficients with broad resonances that appear when 
d-wave channels are included can be observed. 

The precise form of short-range potential and disper- 
sion coefficients are not known. A few percent uncer- 
tainty is not unrealistic. For this Letter we have con- 
structed potentials that lead to a positive _B = scat- 
tering length a for "'^^'''Dy atoms in the m = —8 state. 
This choice is suggested by the recent observation of a 
Bose condensed gas of ^^^Dy atoms at nearly zero mag- 
netic field. It thus possesses a positive scattering length 
at this field (8|. Moreover, we chose the scattering length 
to be approximately equal to the mean scattering length 
P5] for (fictitious) scattering of a van der Waals poten- 
tial with a Cg coefficient equal to the isotropic dispersion 
coefficient. 
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FIG. 4: Scattering length as a function of magnetic field for 
the bosonic isotopes ^^^Dy (top panel) and ^''"Dy (bottom 
panel) . The magnetic state and number of included (. is as in 

Fig-i 



To further elucidate the effect of anisotropy, Fig. [3] 
shows the scattering length of 771 — —8 ^^"'Dy collisions as 
a function of magnetic field when parts of the anisotropy 
are turned off. The top panel displays the case when 
all interactions are included corresponding to the elastic 
scattering described in Fig. [2] The bottom two panels 
show the effect of turning of the anisotropy in the disper- 
sion and magnetic dipole-dipole interaction, respectively. 
The resonance spectra in the three panels are quite dis- 
tinct. The number of resonances differs and, with one 
exception, the resonances are narrower. 

Finally, Fig. |4] illustrates the effect of changing to dif- 
ferent bosonic Dy isotopes. Since to good approximation 
Born-Oppenheimer potentials do not depend on the iso- 
tope, we have solved the coupled-channels equations us- 
ing the appropriate reduced mass /x^- This observation 
has been used in understanding relationships between 
scattering lengths of isotopic combinations of spin-less 
ytterbium [29] , while its limitations for Lithium Feshbach 
resonances have been studied in Ref. [3D]. The field de- 
pendence of the scattering length changes from ^^'^Dy to 
the ^^^Dy isotope. Measurement of resonance locations 
in different isotopes will be invaluable in understanding 
the scattering of dysprosium. 

Conclusion. Applying a full coupled-channels calcula- 
tion for ultracold atom-atom collisions, we have shown 
that the origin of Feshbach resonances in interactions 
between ultracold rare-earth atoms with large magnetic 
moments result from strong scattering anisotropies. Con- 
sequently, by tuning an applied magnetic field we pre- 
dict that it will be possible to observe resonances and 
control collisions even for atoms with zero nuclear spin. 
This study is the first predictions of a Feshbach resonance 
spectrum for rare-earth atoms. 

We have investigated the effects of different short-range 
and long-range anisotropic potentials as well as different 



isotopes on the scattering length of ultracold Dy atoms 
as a function of magnetic field strength. To optimize the 
potentials we must await experimental observations of 
resonances from multiple isotopic combinations. 
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